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Are we missing baryons in Galaxy Clusters? 
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ABSTRACT 

The recent constraints on the cosmological parameters put from the observations of the 
WMAP satellite limit the cosmic baryon fraction in a range that is larger than, and marginally 
consistent with, what is measured in galaxy clusters. This rises the question whether or not we 
are considering all the ingredients of cluster baryonic budget. Carefully weighing the baryons 
in X-ray emitting plasma and stars in cluster galaxies, I conclude that the cluster baryonic pie 
is made by 13 (with a la range of 8-19) per cent of stars, 70 (56-89) per cent of intracluster 
hot medium and 17 (0-33) per cent, and a probability of 73 per cent of being larger than 0, of 
"other" baryons, presumably in the form of warm (10 5 — 10 7 Kelvin) material. 

Key words: galaxies: cluster: general - galaxies: fundamental parameters - intergalactic 
medium - X-ray: galaxies - cosmology: observations - dark matter. 
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1 INTRODUCTION 

The recent analysis of the angular power spectrum of the Cosmic 
Microwave Background (CMB) obtained from WMAP (Bennett et 
al. 2003) has provided constraints on the cosmological parameters 
(Spergel et al. 2003) that confirms with greater accuracy the current 
energy density of the Universe to be comprised by about 73 per 
cent of dark energy and 27 per cent of matter, mostly non-baryonic 
and dark. In particular, the quoted constraint on the baryon den- 
sity, fib, is 0.0224 ± 0.0009/iJ~ 2 , and on the total matter density, 
fi m , is 0.135to QQ9/iJ" Q. Consequently, the cosmic baryon fraction, 
fib/ft m , is equal to 0.166± ;^3 and the ratio between baryon and 
cold dark matter density, fi c = f2 m — fib, is equal to 0.199lQ Qjg. 
These values are expected to be maintained in regions at high over- 
densities that collapse to form galaxy clusters. 

The clusters baryon budget is composed mainly from the X- 
ray luminous baryons, M gas , of the intracluster medium (ICM) that 
becomes hotter upon falling into the cluster dark matter halo by 
gravitational collapse. Other contributions come from the baryonic 
stellar mass in galaxies, A/ ga i, and from other "exotic" sources, 
like intergalactic stars and a still poorly defined baryonic dark mat- 
ter. Given the large uncertainties on the relative contribution from 
baryons that are not accounted for in either A/ gas or M ga i, I qual- 
ify these as "other baryons", M b, as already done in a previous 
work (Ettori 2001) in which I discussed the constraints on clus- 
ter baryon budget from BOOMERANG and MAXIMA-I data. The 
tighter constraints on the cosmological parameters provided from 
WMAP allow now more firm conclusions. 

Therefore, one can put the following relation between the rel- 
ative amount of baryons in the Universe and in clusters with total 
gravitating mass, M to t: 



where / g . 



M gaa /M tot , / gaJ 



M ga .i/Mtot (P 



Ob 



Mb 
Mto 



YBC 



+ 



B 



+ 



fob 

~B' 



(1) 



0.010l°;°°;j h'ol in Fukugita et al. 1998), f oh = M ob /M tot , and 
Y is the parameter representing the cosmic depletion of baryons at 
the virial radius with respect to the global value (~ 0.92 ± 0.06 
from the hydrodynamical simulations of the Santa Barbara Project; 
Frenk et al. 1999). I parametrize the uncertainties on the measure- 
ments of the total gravitating mass and gas mass through the factors 
B and C, respectively. These factors act to increase the total mass 
estimates (i.e. B > 1) if corrections to the hydrostatic equilibrium 
equation are required for bulk motions of the ICM or non-thermal 
pressure support, and to lower the true gas mass (i.e. C > 1) if 
dumpiness is present in the ICM that is assumed to be smoothly 
distributed (e.g. Mathiesen et al. 1999). 

In this Letter, I will analyze the equation Q(i) to assess the 
consistency between the cosmic and the cluster baryon budget and 
(ii) to put significant constraints on / Q b . In a consistent way, I adopt 
the WMAP results on the Hubble constant, H , of 711°,;°^ km s" 1 
Mpc -1 to rescale all the measured quantities. 



2 THE COSMIC BARYON BUDGET FROM WMAP 

Figure Q shows the allowed la region from the WMAP re- 
sults on the cosmic baryon budget with respect to the observed 
gas and baryon fraction for a sample of relaxed galaxy clus- 
ters both at low and high redshift. All these values are esti- 
mated at the overdensity of 200 with respect to the critical den- 
sity (i.e. within the cluster region that numerical simulations 
show to be virialized in a fi m independent way; see, e.g., Evrard 
et al. 2002). Using a Bayesian approach (Press 1996), I mea- 
sure the average gas (baryon) fraction and confirm that is con- 
sistent between the two samples: 0.107l^ o g(0.133±^ o ^) and 
0.111±g;g|| (0.133±g;g|^) (the error-weighted means of the gas 
fraction distribution are 0.116 ± 0.005 and 0.111 ± 0.010 in the 
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Figure 1. Gas (dots) and baryon (diamonds) fractions estimated at an overdensity of 200 for a cosmology fixed to the best-fit WMAP results (Spergel et 
al. 2003). These results constrain the cosmic baryon budget to the shaded region. The central (Her) value of the observed distribution is obtained with a 
Bayesian method (see Press 1996) and is indicated with a dashed (dotted) line. The thickest dashed line indicates the estimate of the baryon fraction corrected 
by the depeletion parameter Y = 0.92. (Left) Low (z < 0.1) redshift sample of galaxy clusters observed with BeppoSAX. The best-fit spectral results were 
deprojected to recover the gas density and temperature profile and to infer the gas and gravitating mass (see Ettori, De Grandi, Molendi 2002). The / ga i 
is estimated from B-band luminosities in Girardi et al. (2002), assuming Af ga i/Ls = 4.5 ± l.Ohioo Mq/L@ (Fukugita et al. 1998). Where Lb is not 
available, the median value of measured / ga i = 0.022 obtained from six nearby systems is used [Lin et al. (2003, Fig. 7) measure / ga j between 0.01 and 0.02 
from K-band luminosity]. (Right) High (z > 0.7) redshift sample of objects observed with Chandra (see Ettori, Tozzi, Rosati 2003). 



B=1.00; C=1.00 B=1.15; C=1.00 B=1.15; C=1.15 




Figure 2. Constraints as function of the depletion parameter Y on the gas (dashed line), stellar (dotted line) and "other baryons" (solid line and probabil- 
ity contours at 1, 2, 3 sigma confidence level) mass fraction normalized to the cosmic value. The shaded region indicates the range of Y permitted from 
hydrodynamical simulations (Frenketal. 1999). 



low— z and high— z sample, respectively. Note that these values 
do not include any correction by the baryonic depletion that is ex- 
pected to raise the gas/baryon fraction by about 8/6 per cent at this 
overdensity. See also Fig.Qand text that follows). These estimates 
are consistent with other, independent, recent determinations (e.g. 
Allen et al. 2001, Pratt & Arnaud 2002) and consistently lower than 
the baryon budget required from WMAP results. It is worth notic- 
ing that only the highest estimates of the gas (baryon) mass fraction 
(e.g. A426, A2142, RXJ1350) are perfectly consistent with WMAP 
results, whereas the other clusters are systematically below them. 



To investigate the systematics that could affect this estimate, 
I change the values of the factors B and C and study their influ- 
ence on the baryon fraction. The factor B, which parametrizes the 
uncertainties on Mtot is expected to be between 1 and 1.15 from 
the cluster mass profiles recovered from both X-ray and lensing 
data (e.g. Allen et al. 2001). The factor C represents the level of 
dumpiness that affects the estimate of M gas in X-ray analysis and 
that simulations show to be lower than 1.2 (Mathiesen et al. 1999). 

Figure|2|shows that higher values of B and C require a more 
relevant role to be played by / Q b, giving a 2a positive detection for 
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Figure 4. Cluster baryonic pie. The la range from the Bayesian calcula- 
tions applied to the eight nearby clusters is shown in parenthesis. A deple- 
tion factor Y = 0.92 ± 0.06 is assumed. The factors B and C are fixed to 
1. We expect larger values of the warm ingredient if either B or C are > 1. 
It is worth noticing that the hot component, / gas , scales as h~ 1,6 and the 
cold component, / ga i, is here independent from h. 

typical value of Y when B and C are 15 per cent larger than the 
null hypothesis of reliable estimates of both Af gas and Af to t from 
X-ray analysis. 

Moreover, I can estimate from the observables the ratio 

i?gas = (n b /fi m - /gal /B) //gas X (Y BC) » (/ob// gaS + 1) 

and to evaluate the probability that R ga , s > 1 and a no-zero value 
of fob is required from the data (see Fig.[3J. 1 obtain values of i? gas 
between 1.3 (B = 1, C = 1) and 1.7 (B = 1.15, C = 1.15), 
with an interval accepted at the 95 per cent confidence level of 0.5- 
2.7. More significantly, this ratio has to be larger than 1 at 76.6 
(B = 1, C = 1) and 92.5 (B = 1.15, C = 1.15) per cent con- 
fidence level. This result gives a high confidence to the conclusion 
that a significant amount of baryons has to be present apart from 
what is observed both at X-ray and optical wavelength. 



3 CONCLUSIONS 

By comparing the recent cosmological constraints from measure- 
ments of the angular power spectrum of the temperature anisotropy 
in the CMB done with WMAP with the observed distribution of the 
gas mass fraction in clusters of galaxies, I conclude that 

(i) galaxy clusters with the highest observed gas/baryon frac- 
tion are well in agreement with the WMAP estimate of the cosmic 
baryon budget. On average, however, a disagreement in the order of 
15-20 per cent is present, the observed cluster baryon fraction be- 
ing a lower estimate of the cosmic one. This implies that estimates 
of the cold dark matter density, Q c , done by applying mean results 
from large sample of objects tend to overestimate it. In this perspec- 
tive, a more "realistic" result is provided from the highest estimate 
in the distribution of the cluster gas (baryon) mass fraction; 

(ii) this dark baryonic component appears to be between a frac- 
tion, and up to 2 times, the measured gas fraction. Values of 
# g as = (/wmap - / ga i)// gas larger than 1 are required from the 
data with a level of confidence of about 80 per cent and more, if we 



are underestimating (overestimating) the total (gas) mass. The most 
probable values are / b //gas = Rg&s — 1 = 0.3 (B = C = 1) and 
0.7 (B = C = 1.15), and lower than 1.7 at 95 per cent confidence 
level. 

It is very unlikely that galactic objects, such as "colored" (red, 
brown, white, beige) dwarfs, stellar remnants and other species of 
MACHOs (see review in Gilmore 1999 and Evans 2003), or in- 
tergalactic ones formed from tidal disruption of cluster dwarfs, like 
planetary nebulae (e.g. Ciardullo et al. 2002), red-giant-branch stars 
(e.g. Ferguson et al. 1998) and supernovae (Gal- Yam et al. 2003), 
can be responsible for such amount of baryons. It is reasonable to 
believe that they can contribute by about 0.2 times / ga i, or 0.04/ gas . 

The most plausible suspect to give so large contribution is 
then a X-ray warm (10 5 K < T < 10 7 K) intracluster medium 
(W-ICM). Large scale cosmological and hydrodynamical simula- 
tions by Cen & Ostriker (1999) and Dave et al. (2001) show that 
the mass fraction at redshift is largely dominated from a warm 
medium, with a relative contribution in mass that is about 2 times 
larger than the amount of hot (T > 10 r K) baryons. However, less 
than 30 per cent of it falls in overdensities > 60 that are typical for 
bound structures in a ACDM universe. Furthermore, Bonamente et 
al. (2002) present evidence of excess in the soft X-ray emission 
between 0.2-0.4 keV in 50 per cent of the 38 clusters in their sam- 
ple of Rosat PSPC observations. They list several suggestions on 
how to explain this excess, originally observed in extreme ultravio- 
let (Lieu et al. 1996), both as thermal and non-thermal component. 
If we assign this emission to the baryons that are lacking in our 
budget, we interpret it as thermal emission due either to the dif- 
fuse/halo component of unresolved X-ray faint cluster galaxies or 
W-ICM. In the first case, we are forced to consider an inexplicable 
large amount of X-ray emitting member galaxies. More plausible is 
then the hypothesis of W-ICM, even though its cooling time tends 
to be very short with the bulk of the radiation in emission line if this 
gas is not primordial. Fabian (1997) suggested that it can be located 
in turbulent mixing layers lying between embedded cold clouds and 
the ICM. However, the traditional picture on the efficiency of cool- 
ing processes in the cluster cores is not supported anymore after 
that XMM and Chandra observations did not report evidence of gas 
cooler than 1-2 keV (e.g. Peterson et al. 2003) and showed a strong 
interplay between ICM, the central active galaxy (e.g. Fabian 2002) 
and merging cool clumps (e.g. Markevitch et al. 2000, Mazzotta et 
al. 2003). On the other side, the production of thermal energy per 
particle due to supernovae related to the star formation activity is 
in the order of 0.4(j)/0.1)(A i 'sNii/10 9 )(10 13 M Q /Mgas) keV for a 
given efficiency r\ in converting the kinetic energy of the explosion 
into thermal energy through galactic winds, and adopting typical 
values of the cluster gas mass and number of type II supernovae as 
required from the observed ICM metallicity. While this energy per 
gas particle is not enough to stop cooling the hot ICM, it can easily 
accommodate for the survival of the warm component. 

Intriguingly, the stronger soft excess detected in Rosat data 
from Bonamente et al. (2002) is measured in objects like A85 and 
A1795, that are the ones lacking most of the baryons with respect to 
the cosmic budget as plotted in Fig.Q Significant detection is also 
present in A2029, A2199 and A3571, whereas a marginal detection 
is associated to A2142. Nevalainen et al. (2003; see also Kaastra et 
al. 2003) confirm with XMM the excess in the soft X-ray emission 
in A1795 and that this excess is best fitted by a thermal compo- 
nent with a characteristic temperature of 0.8 keV, which is about 
an order of magnitude higher than what required from Rosat data 
but still consistent with our energetic arguments. Using their esti- 
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mation of the atom density of the W-ICM in the core of A 1795, 
and assuming that it is broadly distributed like the ICM, one can 
infer a /w-lCM//gas = ^gas — 1 ~ 0.43. In general, values of 
/w-iCM//gas between 0.1 and 0.5 are expected. 

To summarize, clusters seem to have similar behavior in ac- 
cumulating the same relative amount of baryons. It is then their 
peculiar thermal history due to the interplay of merging actions 
and/or activity of the central active galaxy that provides the bary- 
onic ingredients and cook the baryonic pie that we taste and show 
in Figure [4] To prepare it, I have considered only the eight nearby 
clusters which provide a more reliable estimate of / ga i and are less 
affected from systematics in the determination of / gas (see discus- 
sion in Ettori et al. 2003). I have also corrected the gas fraction by 
the depletion factor Y ~ 0.92. The baryonic pie is then made of 70 
per cent of hot ICM, with la range between 56 and 89 per cent and 
a distribution of the calculated / gas //wMAP that spans between 
28 and 143 per cent at 2a confidence level (higher upper limit ob- 
served in A426 that has a most probable / gas //wMAP value of 1 15 
per cent). The cold, stellar component is responsible for 13 (la: 8- 
19) per cent, with an observed distribution in the sample between 2 
and 37 per cent (2a lower and upper limit, with the latter reached 
in A2199, which has a central value of 21 per cent). Finally, a third 
ingredient, probably a warm ICM, contributes by about 17 (0-33) 
per cent (and a probability to be larger than of 73 per cent) with 
a distribution that goes from -29±15 per cent in A426 to 40±12 
per cent in A1795, one of the objects with the largest detected soft 
excess (Bonamente et al. 2002, Kaastra et al. 2003). 
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